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ABSTRACT 

The Weibull modulus, fracture toughness, and thermal properties of a 
silicon nitride material used to make a gas turbine combustor were 
experimentally measured. The location and nature of failure origins resulting 
from bend tests were determined with fractographic analysis. The measured 
Weibull parameters were used along with thermal and stress analysis to 
determine failure probabilities of the combustor with the CARES design code. 
The effects of data censoring, FEM mesh refinement, and fracture criterion 
were considered in the analysis. 

INTRODUCTION 

The need for more fuel efficient transportation and lower engine 
emissions has made ceramics increasingly important structural materials. 
Ceramics offer many beneficial properties, such as high-temperature strength, 
low density, and low thermal conductivity. 

Ceramics, however, are brittle materials. Weibull statistics (1) are 
commonly used for reliability analysis of components fabricated from such 
materials. The calculated failure probability is dependent on the stress 
state and the material properties of the component. It is well-known, 
however, that this procedure is inadequate to describe experimental data 
obtained in multiaxial stress states (2). 

The application of ceramics as structural components requires a more 
rigorous reliability theory. Such a theory should combine linear elastic 


fracture mechanics with established statistical failure theories. Several 
applicable theories have been incorporated into the public domain, Ceramics 
Analysis and Reliability Evaluation of Structures (CARES) ( 3 ) code developed 
at NASA Lewis. This post-processor code, when combined with a finite element 
method (FEM) stress analysis, calculates fast fracture probability of a 
brittle, monolithic structural component. CARES also calculates Weibull 
statistical parameters using data from simple mechanical tests such as 
four-point bending or tensile strength tests. 

The objective of this paper is to demonstrate an analytical method for 
determining the mechanical properties and reliability of a ceramic component. 
The example component chosen for this analysis was a silicon nitride combustor 
designed for use in a ceramic gas turbine (CGT) engine. 

COMBUSTOR PROCESSING AND MATERIAL PROPERTIES 
Material Processing 

The silicon nitride material used to manufacture the combustor was 
Kyocera SN 251 . The combustor was processed by the general flow diagram shown 
in Fig. 1. A mixture of silicon nitride (Si 3N4 ) powder, sintering aids and a 
polymeric binder were milled. The resulting slurry was dried, granulated and 
cold isostatical'ly pressed to form a cylindrical tube. The tube was green 
machined to near final shape, dewaxed, and sintered in a nitrogen atmosphere. 
Final machining produced the combustor shown in Fig. 2. 

Mechanical Properties 

Mechanical properties necessary for the FEM and reliability analyses were 
experimentally determined. The parameters required for stress analysis 
included Young's modulus, Poisson's ratio, thermal conductivity, thermal 
expansion coefficient and specific heat. The CARES analysis requires strength 
data (from which the Weibull parameters are calculated) and Poisson's ratio. 
Fracture toughness was also measured. Although this parameter is not required 
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for predictions made by CARES, it is used in life-time prediction and is 
related to the properties used in the reliability analysis. Strength 
measurements, which are the basis of reliability predictions, are controlled 
by both the flaws inherent in ceramic materials and the fracture toughness. 
Toughness represents the ability of a material to tolerate flaws and thus can 
be considered in long term applications. 

Young's modulus and Poisson's ratio between room temperature (RT) and 
1400 °C were determined by a resonance method in accordance with ASTM C848-78 
(4). For comparison, room temperature Young's modulus was also determined by 
strain gaging a four-point bend specimen. Poisson's ratio was also determined 
with an ultrasonic pulse method at room temperature. Thermal expansion was 
measured as a function of temperature via a differential transformer. Thermal 
conductivity and specific heat were obtained as a function of temperature via 
the laser flash method (5). 

Test specimens for strength and fracture toughness measurements were 
machined from the combustor liner and heat treated in air at 1150 °C for 
1.5 hr to eliminate machining damage. The flexure strength was determined in 
four-point bending in accordance with JIS (6) at temperatures of 25, 1000 and 
1371 °C. The specimens measured 3 x 4 mm in height and width and the inner 
and outer spans were 10 and 30 mm, respectively. A minimum of 27 strength 
tests were run at each temperature. Flaw type and location (connected to the 
surface or completely within the volume) were determined with optical 
microscopy or scanning electron microscopy as necessary. 

Fracture toughness was determined with the single-edge-precracked-beam 
(SEPB) (7) and chevron-notch (7,8) methods in three-point bending. Specimens 
were 3 x 4 x 30 mm in thickness, width and span. The SEPB specimens were 
tested at a stroke rate of 0.5 mm/min. Chevron-notch specimens were tested at 
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0.005 mra/min. The low stroke rate used for chevron-notch testing was required 
to insure stable crack extension. 

Results and Discussion 

Young's modulus, Poisson's ratio, thermal expansion, thermal conductivity 
and specific heat are illustrated in Figs. 3 to 5 as functions of temperature. 
Young's modulus exhibited a continuous drop with temperature. The room 
temperature, mechanically and sonically determined values were within 3 percent. 
Poisson's ratio was almost constant. The thermal expansion increased linearly 
with temperature, resulting in a coefficient of 3.1xlO _6 /°C. The thermal 
conductivity exhibited a gradual drop with increasing temperature, while the 
specific heat increased gradually with temperature. 

A summary of bend strength test results is given in Table I. Both the 
characteristic strength, a Q , and Weibull modulus, m, decreased with 
temperature. The characteristic strength is the stress at which 63.2 percent 
of the specimens failed. This strength parameter is dependent on specimen 
volume and should not be confused with the scale parameter which is 
independent of volume. 

Failure originated from either a large, 30 to 100 pm grain with a 
hexagonal cross section, or a 15 to 70 pm diameter pore. Examples of failure 
origins are shown in Figs. 6(a) and (b). 

The population of volume type failures observed at 25 and 1000 °C was 
insufficient for statistical characterization. However, nearly 50 percent of 
the specimens failed from volume flaws at 1371 °C. This change from surface 
to volume failure indicated healing of surface connected flaws. Six specimens 
out of 31 tested at 25 °C and three specimens out of 27 tested at 1000 °C 
failed from processing flaws connected to the beveled edges of the specimens. 
These failures were classified as surface failures. 
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Fracture toughness results are summarized in Table II. The chevron-notch 
specimens exhibited nonlinear load-displacement diagrams, indicating stable 
crack extension. Fracture toughness, calculated from maximum load and minimum 
stress intensity coefficient, decreased with temperature up to 1200 °C. 
However, at 1371 °C the load-displacement diagrams became severely nonlinear 
and the toughness appeared to increase substantially. Evidently, the 
combination of slow stroke rate and high temperature resulted in creep 
deformation instead of stable fracture. 

Room temperature fracture toughness measured with the SEPB was comparable 
to that of the chevron-notch. However, SEPB measurements at 1371 °C were much 
lower and comparable to the value measured with the chevron-notch at 1200 °C. 
The SEPB produced fast fracture at 1371 °C, while the chevron-notch allowed 
creep deformation instead of stable fracture. 

COMBUSTOR ANALYSIS 
FEM Stress Analysis 

Before performing the reliability evaluation with CARES, both the 
mechanical and thermal stresses in the CGT combustor were analyzed. The 
thermal stresses were especially significant for the combustor. Empi rical ly , 
the critical stress of the combustor was expected to occur around the dilution 
port due to cooling. When an engine is shut-down, fuel is cut off and the 
combustor cools rapidly because of the injection of fresh air. The resulting 
thermal shock phenomenon produces high, transient tensile stresses on the 
surface. The FEM analysis was designed to simulate the worst transient 
thermal stress in an emergency stop. 

Figure 7 shows the finite element model of the combustor. Exploiting 
symmetry, only 1/16 of the combustor was discretized. Table III shows gas 
temperatures and heat transfer coefficients for both steady state and 
transient conditions. 
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temperatures and heat transfer coefficients for both steady state and 
transient conditions. 

Gas temperature during running was set at 1350 °C for the inside of the 
combustor and at 913 °C for the outside and the primary hole. The heat 
transfer coefficient was assumed to be the same for the steady state and 
transient conditions. 

Figure 8 shows the steady state temperature distribution obtained by the 
FEM analysis. There was less than a 70 °C temperature difference in the 
combustor. Thus, the resulting stress level was relatively low, with a 
maximum of 60 MPa located at the edge of the primary hole. 

During the transient stage, however, the maximum principal stress of this 
part changed with time and reached a peak value of 163 MPa after 30 sec, as 
shown in Fig. 9. Figures 10 and 11 show the temperature distribution and the 
maximum principal stress, respectively, at 30 sec. There was almost a 200 °C 
temperature difference between the primary hole and the bottom of the 
combustor. This large temperature difference resulted in the high tensile 
stress at the edge of the primary hole. 

These results indicated relatively low stresses at steady state and 
higher stresses after combustor shut down. Thus, the reliability analysis in 
this study concentrated on the tensile stress state at 30 sec after shut down. 
Description of the CARES Program 

The CARES program is based on the theory proposed by Batdorf (9) in which 
linear elastic fracture mechanics (LEFM) is used to derive multiaxial effective 
stress equations. The critical stress of a crack is defined as a remote, 
uniform uniaxial stress which causes fracture when applied normal to the crack 
plane. With the assumption that the cracks are randomly oriented, the 
effective stress equations are used to calculate the probability that a crack 
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with a given critical stress level is oriented with respect to the applied 
loads such that it would cause fracture. This orientation probability is 
multiplied by a crack existence probability, namely the probability that a 
crack with a given critical stress exists. 

For structural applications, integration must be performed over the 
volume of the component to account for the variation in stress with location. 
In CARES, the results of finite element analysis are used for this purpose. 

To apply the Batdorf theory to a practical problem, one must make 
simplifying assumptions about the type of flaws present in the material. In 
the CARES program, several different flaws types are available. It is usually 
recommended that penny-shaped and semicircular cracks be used in the volume 
and surface analyses, respectively, as they are more representative of small 
embedded flaws. 

Further assumptions must be made about the condition that will lead to 
fracture. The LEFM criteria available in CARES include the normal stress 
criterion, the total strain energy release rate criterion, the maximum 
principal stress criterion, the principle of independent action ( P I A ) 
criterion, and a mixed-mode interaction equation called Shetty’s criterion 
(10). Shetty's criterion can be used with a constant, C, defined as the ratio 
of fracture toughness in pure mode II loading and the mode I fracture 
toughness. The constant can be selected to fit data from mixed-mode fracture 
toughness testing, or selected so that this criterion simulates out-of-plane 
crack extension theories, including the maximum strain energy release rate 
(Gmax) criterion. 

The material parameters which describe the variability in strength must 
be obtained from fracture strength measurements. The measurements usually 
have a mixture of volume and surface failure origins. The locations of these 


7 


origins may be distinguished (censored data) or ignored (uncensored data). The 
distribution of strength is often described with Weibull modulus, m, and 
characteristic strength, a Q. CARES can estimate these parameters by either 
maximum likelihood (ML) or a least squares fitting (LS) technique. 

For the reliability analysis, the characteristic strength is not directly 
used, as it depends on the effective volume or surface area of the test 
specimen. Instead, the value corresponding for a uniformly stressed specimen 
of unit volume or area, the scale parameter, is used. The scale parameter is 
theoretically independent of the material volume tested, however, it must be 
calculated from the characteristic strength and Weibull modulus and depends 
upon the assumption that the Weibull weakest link theory is valid. 

For the Batdorf theory, the probability of a crack existing is assumed to 
be proportional to a crack density function coefficient, kg. This coefficient 
is determined from the characteristic strength by equating the Batdorf and the 
Weibull failure probabilities for the given test conditions. It should be 
noted that the value of kg depends on the fracture criterion that is used in 
this calculation. There are two methods available in CARES for calculating 
kg. The first one assumes that in uniaxial loading only mode I fracture can 
occur. This restriction is due to the assumption that a sufficient number of 
large cracks are oriented perpendicular to the loading direction, and that 
these will preferentially fail. This method results in a shear-insensitive 
crack density coefficient. The second technique assumes that the selected 
mixed-mode fracture criterion is equally applicable for uniaxial and 
multi axial loading. This shear sensitive method will always lead to a lower 
crack density coefficient and a lower resultant multiaxial failure probability. 
Failure Probabilities 

Statistical parameters . The rationale to be used in the selection of 
material parameters for reliability calculations is not straight forward. 
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Goodness-of-f i t tests provide some indication of whether the ML or LS 
parameters are better suited for the data at hand. 

Our rationale used in choosing material parameters for the reliability 
calculations was to provide conservative values of failure probabilities of 
the component. In the low failure probability regime, a lower value of the 
Weibull modulus results in higher failure probabilities. Failure probabilities 
are not strongly dependent on the characteristic strength, whereas the Weibull 
modulus has a much more significant effect. In accordance with the 
aforementioned arguments, LS parameters were chosen for the analyses. 

In this case, the lack of volume flaw data at RT and 1000 °C caused both 
ML and LS parameters to be questionable, resulting in undependable 
goodness-of-f i t tests. Thus, the volume flaw characteristic strength and 
Weibull modulus at 1000 °C were obtained by graphical interpolations of the 
corresponding values at RT and 1371 °C, as shown in Fig. 12 and listed in 
Table I. 

For comparison, the data was also analyzed without censoring, i.e., 
fracture origins were not considered. Obtained parameters are listed in 
Table IV. 

Results and discussion . The failure probabilities resulting for several 
fracture criteria are listed in Table V. Very small failure probabilities were 
obtained for all criteria. Failure probabilities were larger for the surface 
analysis than the volume analysis in all cases. The maximum total failure 
probability was 2.463xl0~ 9 for Shetty's semi-empirical case with C = 0.80, and 
the minimum failure probability was 6.962xl0 -19 for the PIA model. 

Parameters obtained from the uncensored data were also used to calculate 
the reliability of the component. For uncensored data, the estimated failure 
probabilities for Shetty's semi-empirical case were I.l06xl0 -13 and 2.315xl0 -12 
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for volume and surface analyses, respectively (see Table VI), using LS 
parameters. Comparing these values with the corresponding ones in Table V, it 
is evident that there are large differences in the failure probabilities for 
censored and uncensored data. Although parameters determined from uncensored 
data may have reasonable statistical characteristics because of the larger 
sample size, differences can arise with such simplifying assumptions. 

Additional calculations of the worst case situation were made by 
selecting the 90 percent lower confidence limits of the Weibull moduli and 
characteristic strengths obtained by ML. The 90 percent limits are not 
available from the CARES program for sample sizes smaller than five. For 
these cases the 90 percent lower bound was estimated using the same factors 
that would have been used for a sample size of five. The failure probability 
increased to 5.302x10“® for this extreme condition, as shown in Table VI. 

The calculated stresses may depend on the FEM mesh used. In this case, the 
calculated stresses are concentrated on the edge of the hole. The differences 
in stress states may be enhanced by raising them to the power of the Weibull 
modulus. Failure probabilities for a Standard (coarse) and a fine mesh are 
listed in Tables V and VI for Shetty's failure criterion. The failure 
probabilities obtained with the fine mesh were smaller than those of the 
standard mesh. 

All of the above calculations resulted in very small failure probabilities 
of the combustor, even for the worst case calculation. Thus, it is concluded 
that this combustor can safely operate under the assumed loading conditions. 
CONCLUSIONS 

Calculations of statistical parameters and failure probabilities for a 
SN251 combustor were performed with the CARES computer code. The failure 
probabilities under extreme operating conditions were presented assuming a 
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worst case scenario. From these calculations it was concluded that this 
combustor can be safely operated for the short-term under the assumed 
conditions. It should be emphasized that a quantitative estimate of the 
failure probability gives good guidance for designing ceramic components. 

The importance of statistical analysis, and censoring of volume and surface 
flaws were demonstrated. Simple bend tests did not provide sufficient volume 
fracture data for this analysis. To obtain more volume fracture data, tensile 
testing may be required. 

In this paper, only the fast fracture failure mode was considered. For 
the lifetime evaluation of such components, additional fatigue testing as well 
as further theoretical work are required. 
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TABLE I . - BEND STRENGTH TEST RESULTS 



Least squares fit 

Maximum 

1 ikel ihood 

Temper- 

Flaw 3 

Number of 

a 9’ 

Weibul 1 

V 

Weibul 1 

ature, 

location 

specimens 

MPa 

modulus 

MPa 

modulus 

°C 







25 

S 

26 

694 

18.36 

694 

18.65 


V 

3 

752 

21.44 

749 

32.44 


u 

2 





1000 

s 

26 

591 

13.78 

592 

13.33 


V 

1 

b 623 

b 17.00 

b 623 

b 23 .80 


U 

0 





1371 

s 

29 

565 

10.12 

574 

8.86 


V 

21 

585 

13.95 

585 

15.31 


u 

0 






a S = surface flaw; V = volume flaw; U = unknown origin. 
^Interpolated values (see section 3 and 3(a)). 


TABLE II. - FRACTURE TOUGHNESS DATA 


Test 

Temperature , 

Number 

Kic, 

method 

°C 

of tests 

MPa'm- 5 

Chevron 

25 

5 

7.9 

notch 

800 

5 

7.1 


1000 

3 

6.9 


1200 

4 

6.0 


1371 

2 

10.4 

SEPB 

25 

7 

7.4 


1371 

4 

6.1 


TABLE III. - GAS TEMPERATURE AND HEAT 
TRANSFER COEFFICIENTS 



Inside 

Outside 

Primary hole 

Heat transfer 
coef f icient 
(W/m 2 K) 

85.5 

113.3 

1032 

Gas temperature, °C 


Steady 

state 

1350 

913 

913 

Transient 

state 

300 

300 

300 


13 


TABLE IV. - UNCENSORED WE I BULL PARAMETERS BY 

THE LEAST SQUARES METHOD 


Temperature, 

°C 

Scale parameter , 
MPax(mm 2 ) 1 

V 

MPa 

Wei bull modulus, 
m 

25 

844 

690 

18.90 

1000 

776 

590 

14.05 

1371 

770 

544 

11.25 


TABLE V. - FAILURE PROBABILITY OF THE COMBUSTOR PREDICTED 

USING THE CARES CODE 


Cr i teria 

Volume 

Surface 

Total 

Shetty's semi-empirical 
Penny shape (semi-circular) 
Mode I fracture 

1.823xl0- n 

2.445xl0~ 9 

2.463x10-9 

Shetty's semi-empirical 
Penny shape (semi-circular) 
Shear sensitive 

3 . 874xl0 -13 

7.719x10-1° 

7.723xl0- 10 

Strain energy release rate 
Penny shape (semi-circular) 
Mode I fracture 

1.146X10- 12 

9.929xl0~ 12 

9 . 940xl0 -10 

Strain energy release rate 
Penny shape (semi-circular) 
Shear sensitive 

4 .217xl0 -13 

7 . 050xl0 -10 

7.055xl0- 10 

Normal stress fracture 

5.223xl0 -13 

7 . OOlxlO -10 

7.006xl0- 10 

Maximum tensile 

Penny shape (semi-ciruclar) 

9.954xl0~ 13 

7 .OOlxlO -10 

7.011x10-1° 

PI A model 

5 .246xl0 -13 

6 . 957xl0~ 10 

6.962x10-1° 
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TABLE VI. - FAILURE PROBABILITIES PREDICTED USING THE CARES CODE 


[Criteria: Shetty's semi-empirical mode I fracture criterion 

penny shape (semi-circular).] 



Volume 

Surface 

Total 

Uncensored 

l.ioexio- 1 ^ 

2 .315xl0 -12 

2 . 426xl0-l 2 

Extreme Wei bull 
modulus 

2.888xl0~ 12 

5. 302xl0- 6 

5.302xl0 16 

Refined mesh 

1.249xlO~ 16 

4.976xl0" 12 

4 . 976xl0 -12 

Refined mesh 
with extreme 

Weibull modulus 

1 .665xl0 -16 

2.482xl0- 8 

2.482xl0- 8 
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Sintering 

additives 


Silicon nitride powder 



Polymer 

binder 



Figure 1 . — Flow diagram for processing the SN251 Figure 2. — View of the SN251 silicon nitride combustor, 

silicon nitride combustor. 
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Figure 3. — Young's modulus and Poisson's ratio 
of SN251 as a function of temperature. 


Figure 4. — Thermal expansion of SN251 as a function 
of temperature. 
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Thermal conductivity, W/m 





Figure 5.— Thermal conductivity and spe- 
cific heat of SN251 as a function of tem- 
perature. 



(a) Large, surface connected grain. 



(b) Volume pore. 


Figure 6. — Failure origins. 
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Figure 7. — Finite element 
mesh of a 1/16 section 
of the combustor. 


Figure 8. — Contour plot of the steady state 
temperature distribution in the combustor. 



Time, sec 


Figure 9. — Maximum principal stress in the com- 
bustor as a function of time. 
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Figure 10. — Temperature distribution (°C) 
at the peak value of the maximum 
principal stress. 


Figure 1 1 . — Maximum principal stress 
distribution (MPa) occurring 30 sec 
after engine shut-down. 
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